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Abstract 

We investigate the motion of a nano tippe top, which is formed from a Cgo fullerene, 
and which is assumed to be spinning on either a graphene sheet or the interior of a 
single-walled carbon nanotube. We assume no specific geometric configuration for the 
top, however for example, the nano tippe top might be formed by joining a fullerene Cgo 
with a small segment of a smaller radius carbon nanotube. We assume that it is spinning 
on a graphene sheet or a carbon nanotube surface only as a means of positioning and 
isolating the device, and the only effect of the graphene or the carbon nanotube surface 
is only through the frictional effect generated at the point of the contact. We employ 
the same basic physical ideas originating from the classical tippe top and find that the 
total retarding force, which comprises both a frictional force and a magnetic force at the 
contact point between the Ceo fullerene and the graphene sheet or the inner surface of 
the single-walled carbon nanotube, induces the Cgo molecule to spin and precess from a 
standing up position to a lying down position. Unlike the classical tippe top, the nanoscale 
tippe top does not flip over since the gravitational effect is not sufficient at the nano scale. 
After the precession, while the molecular top spins about its lying down axis, if we apply 
the opposite retarding magnetic force at the contact point, then the molecule will return 
to its standing up position. The standing up and the lying down conflgurations of the 
nano tippe top during the precession and retraction processes demonstrate its potential 
use as a memory device in nano-computing. 

1 Introduction 

The discovery of fullerenes [l| and carbon nanotubes 0] has led to numerous studies on 
their properties and their various potential applications in nano devices. In this paper, 
we focus on the mechanics of a nanoscale tippe top comprising a fullerene Ceo which 
for example might be joined to a carbon nanotube and is spinning on a graphene sheet 
or inside a large single-walled carbon nanotube. We assume that the nano tippe top is 
at the equilibrium configuration either on the graphene sheet or inside the outer carbon 
nanotube and that the spinning occurs due to the application of an external magnetic 



'Corresponding author's email address: yc321@uow.edu.au 



1 



field. We refer to this structure as a nano tippe top or simply as a nano top. We note 
that in the case of the top spinning inside a nanotube, that although the effect of the 
outer carbon nanotube is not directly incorporated into our calculations, we have in mind 
that the outer carbon nanotube acts as a barrier between the nano top and the external 
environment. 

The flip over of the classical tippe top has attracted much attention due to its "ap- 
parent" violation of the principle of conservation of energy during the top's inversion (i.e. 
the rise in its center of mass results in the sudden increase in its potential energy but 
apparently no gain in other energies) as well as the lack of a complete mathematical de- 
scription of this inversion phenomenon 0, 0, S, 0, 0, @, 0] • However, Cohen [13] provides 
an analysis and a numerical study on the tippe top. He views the tippe top as an eccentric 
sphere, for which its center of mass is different from its geometric center, and incorporates 
into his model Coulomb friction at the contact point to describe the tippe top's motion. 
For example, this eccentric sphere can be manufactured either by creating the top with 
two different mass densities or by puncturing the sphere and inserting a stem. At the 
nano scale, the latter method may be achieved by introducing a defect on the surface of 
a fullerene Cgo and then joining to this defect a small segment of a carbon nanotube. We 
refer the reader to Nasibulin et al. [11] for the possible creation of such nano tippe tops. 
Classically, the friction between the eccentric sphere and the relatively rough surface plays 
a vital role for the top's inversion, and the sudden gain in the potential energy during 
the inversion results from the loss in the top's rotational kinetic energy, which can also 
be observed from the reduction in the spinning of the inverted top. Furthermore, since 
gravity and the normal force act solely along the axis fixed in space and provide zero 
torque to the top, the friction, which acts offset from the space-fixed axis and against the 
rotational motion of the top, and is therefore the only source providing an external torque 
to slow the spinning of the top down 0]. In addition, Ueda et al. 13] show theoretically 
that the top's initial spinning and the ratio of the top's two different principal moments 
of inertia at the center of mass play an important role to determine the top's inversion. 
For example, the flip over phenomenon can only occur when its initial spinning is above 



a certain threshold. Since friction at the nano scale is ultra- low [14l . Il5l ]. our numerical 
results show that the friction itself is insufficient to make the nano tippe top precessing 
about the fixed-body axis. Therefore, in this paper we introduce a sufficiently large re- 
tarding magnetic force at the contact point to act analogously to the effect of friction. 
However, we find numerically that instead of hipping over, the nano top prefers to spin in 
a stable lying down configuration, which suggests that the effect of gravity is negligible at 
the nano scale. In particular, the nano top behaves more like a hard-boiled egg spinning 
on a rough surface [ifii] . 

As mentioned above, to induce the spinning effect of a nano top we introduce a re- 
tarding magnetic force at the contact point between the fullerene Ceo and the carbon 



nanotube's wall. We find from Wood et al. [171 ] who experimentally produce ferromag- 
netic fullerenes Cgo that the magnetically strongest fullerenes are formed at 800 K with 
the magnetic moment per molecule of 0.38^B; where denotes the Bohr magneton con- 
stant. This result indicates that we can initiate the spinning of a fullerene in a preferred 
direction by applying an external magnetic field to the center of the fullerene and similarly 
at the contact point for generating the retarding magnetic force. 

In the following section, we briefly state vector equations which are used to describe 
motion of a tippe top. We note that explicit forms of these equations are given in Appendix 
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A. In Section III, we provide numerical results for a nano tippe top, which is driven by a 
constant magnetic force in both x- and y-directions. We verify our numerical schemes by 
examining the classical tippe top, as presented in Appendix A. Conclusion of the paper is 
given in Section IV. In Appendix B, we use the basic equations as given in Appendix A 
to study the stability for certain configurations of a nano tippe top and in particular 
Appendix C considers the compatibility between our numerical results and asymptotic 
expansions when the nano top is in the lying position. Finally, while Section III considers 
Hx = Hy = H, in Appendix D we assume a magnetic force which is applied only in the 
y-direction and only for a finite time to- 



In this section, we state the equations of motion for a nano tippe top, which is schemat- 
ically illustrated in Fig. [H Owing to the axially symmetry of the density of the nano top 
and its eccentric structure, we assume that the center of mass O is at a distance a away 
from its geometric center C. We denote coordinates {X, Y, Z) as the space-fixed frame 
of the top while {x, y, z) as the body-fixed frame. We adopt (0, 9, -0) as the usual Euler 
angles relative to the space-fixed frame so that $7 = </>. In addition, we choose the axes 
OY and Oy such that they coincide for all time. Therefore, the space-fixed frame and 
the body- fixed frame are different only through the rotation of the nutation angle 9. Fur- 
thermore, we suppose that the spinning about the space-fixed frame is ft{t) = (0,0, $7) 
and the spinning about the body-fixed frame as uj where il. denotes angular frequency 
about the Z-axis, and the nano top is assumed to be initially spinning about the axis Oz 
with an angular frequency uj. The distance between the surface and the center of mass is 
h{9) = R — acos9 and the position vector OP is r = (asin^,0, —h[9)). From the above 
quantities, the total angular velocity n and the total angular momentum L of the nano 
tippe top's system in the body-fixed frame are given respectively by 



2 Equations of motion 



z 

As 




Figure 1 ; Schematic of nano tippe top 




L = {{Clo - An COS 9) sin 9, A9, AQ sin^ 9 + Cu cos 9) 



(1) 
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where {A, A, C) denote the principal moments of inertia at the center of mass O. In 
addition, the translational velocity at the contact point P with respect to O can be 
written as Vrp = n x r and hence the sliding velocity at P is given hy Vp = X + v^p, 
where X = {ux,Uy,Uz) denotes the velocity of the center of mass of the nano top. The 
equations of motion of the nano top can then be determined for the six degrees of freedom 
motion comprising three rotational equations derived from the Euler equations and three 
translational equations derived from Newton second law, namely 

dL 

— +nxL = rx{N + F), 

m{X + n X X) = N + F + W, (2) 

where m, N, F, W denote the mass, the normal force, the retarding force and the weight 
of the nano top, respectively. Explicit forms of these equations of motion are given in 



Appendix A of Ueda et al. [131] and they are also briefly stated in the appendix of this 
paper. Next, we determine a suitable form of the frictional force for the proposed nano 
tippe top system. Various theoretical [l^, 19] and molecular dynamics studies [i^, 21] 



suggest that under the low velocity limit, the frictional force between two molecules is 



linearly proportional to their relative velocity. In addition, Heo et al. |22l ] propose that 
the frictional force between a fullerene and a carbon nanotube is also proportional to the 
fullerene's normal reaction, namely \F f \ = fiN. Further, molecular dynamics simulations 
of Heo et al. [2?] show that the frictional coefficients for various nanostructures, namely 
single-walled carbon nanotubes, nanopeapods and double walled carbon nanotubes, under 
a low pressure regime, are essentially the same, i.e. fi = 0.13. In this paper, we also incor- 
porate a retarding magnetic force to the model and we propose that the total retarding 
force of the nano top inside the carbon nanotube is given by 

F = -Ff-H{B) = -^^^vp-H{B), (3) 

1 '^P 1 

where H{B) denotes the retarding magnetic force acting at P and vp \s the velocity at 
the contact point P. 



3 Numerical results and discussion 

A fourth order Runge-Kutta method [2^ is adopted here to numerically solve this system 
of six ordinary differential equations, namely Eq. ([2]). Ueda et al. [l^] show that in macro 
scale the top's inversion is strongly dependent on the ratio of two principal moments of 
inertia at the center of mass. Therefore, we check whether the inversion of the nano top is 
possible by running the numerical scheme with various values of C/A G (0, 1), but we find 
that this ratio does not effect our numerical results. However, in all cases a sufficiently 
large initial spinning of the nano top is required. Therefore, we choose the following 
physical parameters for the numerical iteration: R = 3.55 A, a = O.li?, m = 1.196 x lO^^'* 
kg, A = {2/?>)mB? , C = 0.5^4 and = Hy = 0.1 zN with the following initial conditions: 
6 = 0.1, = 0, = 0, a; = 100 and X = (0,0,0). These initial conditions can be 
interpreted that we release the nano top with its initially spinning (100 Hz) about the 
z-axis, having 0.1 rad defiection from the Z-axis and zero sliding velocity at the contact 
point P. We use 500 grid points to carry out the numerical iteration and the numerical 
results obtained for 9, Q. and uj are illustrated in Figs. [21 [3] and HI respectively. 
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Figure 2: Nutation angle 9 for g 0, 9.8 and 100 ms during precession 




time(s) ^10-° 

Figure 3: Angular frequency about the Z-axis during precession 

Under the total retarding force, which comprises both the frictional and the magnetic 
forces at the point of contact P, Fig.[2]shows that the nano top precesses from its standing 
up configuration {6 = 0) during the first 5 /is but asymptotically approaches the lying 
down configuration (9 = tt/2) after 50 ^s with the decrease in its oscillating amplitude. 
During the precession, as shown in Fig. [3l the angular frequency about the Z-axis 
increases dramatically in the first micro second indicating the sudden drop down of the 
nano top but oscillating around zero while the nano top has laid down. On the other 
hand, from Fig. S] the angular frequency uj about the z-axis also drops from its initial 
positive value to zero and then monotonically increases in magnitude to 30 MHz but in 
the opposite direction with respect to the direction of the initial spinning. This implies 
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Figure 4: Angular frequency cu about the z-axis during precession 




Figure 5: Precession of nano tippe top 



that the nano top reverses its spinning direction and gains its spinning speed due to the 
opposite motion and the energy gained from the retarding magnetic force respectively. 
We also note that the center of mass of the nano top moves in the X, Y-directions only 
resulting in sliding friction but remains intact in the Z-direction (see Figs.llOland [TT]) . The 
whole precession process is illustrated in Fig. [5l Unlike the classical tippe top, we observe 
no flip over phenomenon in the nano top's precession. This no inversion phenomenon 
arises from the fact that the gravitational force is negligibly small at the nano scale. As 
confirmed by Fig. [21 the precession is not effected by the gravitational force since the 
numerical results for g = and 9.8 ms~^ almost coincide with each other, and the same 
behaviour is obtained for g = 100 ms~^. 

Once the nano top is spinning about the lying down axis, it is important to determine 
a possible way for which the nano top can retract to its initial standing up position. We 
find that the nano top will retract smoothly back to its standing up axis upon applying 
a magnetic force of the same magnitude but in the opposite direction to the previous 
retarding magnetic force and the numerical solution for 9 is shown in Fig. [6l 

Accordingly, by adjusting the magnetic field, we can manually control the nano top's 
switching between the standing up and the lying down states and hence it can be utilized 
as a nano-computing memory device. Moreover, the main advantages of utilizing the nano 
tops as a memory device are that it possesses a remarkably short relaxation time (~ 50 
us) resulting in a higher computational speed, and it is small in size and hence provides a 
larger memory capacitance. Most importantly, it is simpler to control as compared with 
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Figure 6: Nutation angle 6 after applying the reversed retarding magnetic field at P 

the electron-spin quantum memories. Finally, the integration of the self-assembled hybrid 
nanostructure known as nanopeapods and the ideas of nano tippe top developed here may 
lead to practical computing memory devices, which currently require large numbers of bit 
handling. 

4 Conclusion 

A nano tippe top formed from a Cgo fullerene and spinning either on a graphene sheet or 
inside a carbon nanotube is investigated as a possible candidate for a computing memory 
device. The equations of motion for such a nano top are described and we find that while 
the retarding magnetic force makes the nano top precess, it does not flip over as in the 
classical tippe top, but due to the fact that gravity is negligible at the nano scale, it adopts 
a lying down position. In addition, while the nano top is in the lying down position, if 
wc apply the magnetic force which is of the same magnitude but in the opposite direction 
to the previous retarding magnetic force, then the nano top will return to its standing up 
position. Hence, the standing up and lying down configurations of the nano top might be 
considered as two bit states, which gives rise to their potential use as a future memory 
device. 
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A Equations of motion 

Here, we briefly state explicit forms of equations of motion obtained from Eq. ([2]). From 
Eq. the three rotational equations are given by 

Ae = -n{Cuj - An cos e) smO - aN sine- h{e)F^, 

A sin eh = {Cuj - 2AQ. cos e)e + {a - R cos e)Fy, 

Cd; = R sine Fy, (4) 

where h{e) = R — a cose. From Eq. ^2 the three translational equations are of the form 

miix = mVtUy + F^, 
miiy = —mnux + Fy, 

mUz = N — mg, (5) 



where 



Fee 



-fiN{u,-h{e)ej 



jti^ - h{e)e^ + {uy + [R{uj - n cos e) + an] sin ey 

■fiN {uy + [R{uj - n COS e) + an] sin 6*} 



\ux - h{e)ey + {uy + [R{uj - J] COS 61) + an] sineY 
-Hy, (6) 

and Hx and Hy denote the strength of the retarding magnetic force in the x- and y- 
directions respectively. By multiplying sine both sides of Eq. ^2 we obtain 

A sin^ eh = {Cuj - 2An cos e) sin ee + {a - R cos e)Fy sin 6*. (7) 

From Eq. we have Fy sine = Cuj/R, which upon substituting into Eq. d?]) we have 

A sin^ eh + 2 An cos esinee = Cuj sin ee - Cdj cos e + aCuj /R, 

which can be written as 

^d{nsin'^e) aC duj ^d{uj cose) 

Thus, from Eq. ([8]) we obtain 

Ansin^e = aCu:/R-Cojcose + J*, (9) 
where J* is a constant of integration. By multiplying R both sides of Eq. ([9]) we have 

J = AnRsin^e + Cuj{Rcose-a), (10) 



where the constant J is referred to as the Jellett's constant [Ij] and our numerical result 
indicates that this constant is zero throughout the precession process. 



As a benchmark, we use our numerical scheme to show the behaviour of the classical 
macro scale tippe top with zero magnetic force. Here, the values of parameters are taken 
to be = 0.015 m, a = O.IR m, M = 0.015 kg, A = C = (2/b)MR^ and /i = 0.1, with 
initial conditions: 6 = 0.1, uj = 100, = and X = (0,0,0). Thus, from Eqs. §i) and 
([6]) we obtain numerical results for nutation angle 9{t) as illustrated graphically in Fig. [71 
From this figure, we can see that 9 approaches vr implying that the top flips over, which 
is consistent with Ueda et al. ul 




time (s) 



Figure 7: Nutation angle for classical macro scale tippe top 



B Stability analysis 



In this section, we use a simple stability argument to investigate the stability of the 
nano top for = 0, 7r/2 and tt. Without the external magnetic retarding force, i.e. 
Hx = Hy = 0, we consider a small perturbation around the standing up axis, namely 
9 = e ioT small e > 0. Eq. (j3|)i becomes 

e = -aNe/A, (11) 

which implies that the spinning nano top is always stable. If we incorporate the external 
magnetic retarding force, i.e. Hx = Hy = H, which is a constant. For 9 = e, Eq. @i 
becomes 

e = -aNe/A + {R- a)H/A, (12) 
which also implies that the nano top is stable. We note that for the case 9 = e, we assume 

n = o. 

Next, we check the stability of the nano top at 9 = ir. Upon substituting 9 = tt — e, 
Eq. @i becomes 

.._ {C + A)n'^ + aN {R + a)H 



noting that uj = Q when 9 tends to vr. Eq. imphes that the nano top is always 
unstable, and therefore it does not flip over. 

For the lying down position, i.e. 9 = 7r/2 — e, Eq. becomes 



An"^ -aH HR- Cujn - aN 
e = : e 



A 



(14) 



which is stable if AQ"^ - aH > or 17^ > aH/A. 



C Asymptotic expansion for 9 = 7t/2 

In this section, we check the compatibility between our numerical results and the asymp- 
totic expansions for the governing ordinary differential equations given in Eq. ([2]) or Eq. ^ 
when the nano top is lying down and the time t is sufficiently large. Upon substituting 
9 = tt/2 into Eqs. (HDi, (gDa, dSDi and (P2, we obtain 



n = RH/{Cuj) 



Ijj 



-RH/C, 



mUr 



mQuy — H, muy 



-rriQux — H, 



(15) 



respectively. Noting here that we assume Hx = Hy = H. From Eq. (|15p 9. upon integrating 
both sides by t, we have 

= -RHt/C + ci. (16) 

This asymptotic expansion and its corresponding numerical result for u are plotted to- 
gether in Fig. [HI We note that the constants used here are given by i? = 3.55 A, C = 0.5^, 
A = (2/3)mi?2, m = 1.196 x lO'^^ kg and H = 0.1 zN. We fit Eq. ([IHD with the numerical 
solution for uj and obtain ci = 3 x 10^. Upon knowing uj, we can determine Q by utilizing 



Numerical result for to 

-(RH/C)t+3x10° 



2 3 

time (s) 



Eq. ([15]) 1 as 



Figure 8: Asymptotic expansion for u; 
RH 



11 



where this asymptotic equation decays to zero for a sufficiently large t and the corre- 
sponding result with its numerical solution are plotted in Fig. [9l Given the asymptotic 



- Numerical result for 
-RH/(3x10''C-RHt) 



time (s) 



Figure 9: Asymptotic expansion for 

expansions for 0, we can determine the asymptotic expansions for both Ux and Uy, which 
are given by 



mu,. 



mRH 
{ciC-RHtfy ~ ' 



mRH 



{ciC - RHt) " 



H. 



We can solve Eq. (jlSp analytically by introducing v = Ux + iuy, where i 
multiplying i both sides of Eq. ()18l) 9 and combining with Eq. (jlSP i we have 



v + 



RHi 



H 



(ciC - RHt) 



[l + i). 



m 



(18) 
T. By 

(19) 



By multiplying both sides of Eq. (fT9]l by an integrating factor exp (— i log (ciC — RHt)) 
l/(ciC - RHty we obtain 



d 



H 



(1 + 



dt \{ciC - RHty J 
which on integrating we deduce 

V (1 + 1 



m (ciC - RHty 



{ciC - RHt)^-' + C2, 



(ciC - RHty (1 - i) Rm 
where C2 denotes a constant. Thus, the solution of Eq. (I19p is given in the form 



y(^t) = (1 + (Clg - RHt) ^ ^^^^lo^(c,C^RHt) 



(l-i) 



Rm 



(20) 
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To determine the complex constant C2 to fit our numerical data, we assign at t = t*, 
Ux{t*) = ui and Uy{t*) = U2, where ui and U2 are constants. Thus, C2 is given by 



where X = ciC — RHt*. Upon substituting into Eq. ()20p and simplifying we have 



where = log {[ciC — RHt]/ X). Next, we expand the above equation using Euler's 
formula, e*^ = cos (3 + i sin (3, to obtain 



Since v = Ux + iuy , Eq. ([2T]) gives rise to analytical solutions for Ux and Uy , namely 



The solutions (f22]) for Ux and are plotted together with their corresponding numerical 
results in Figs. llOland llll respectively. We note that the values of constants used here are 
taken to be t* = 5 x 10~^ s, ui = —6 x 10"'' ms~^ and U2 = —4 x 10"'^ ms~^. As can be 
seen from these figures, the solutions (|22p agree with the numerical results for sufficiently 
large t. 



D Retarding magnetic force as step function 



In this section, we investigate the precession of the nano top subject to a magnetic force 
which is applied only in the y-direction and only for a finite time to- In particular, we 
consider Hy = H}i{tQ — t), where H is a constant representing the strength of the magnetic 
force, H(t) is the Heaviside unit step function and Iq denotes the time when the magnetic 
force is switched off. Two cases are examined, namely to = 0.8 x 10"^ and 10~^ seconds. 
We observe that in the former case the nano top precesses from its initial standing up 
position to 9 = 2.1 soon after the retarding magnetic field is switched off and then it 
oscillates about 6 = tt/2, as demonstrated in Fig. [T2j The effect of the Heaviside function 
can be seen from the behaviour of a;, which is shown in Fig. [131 From this figure, before the 
switch off time to = 0.8 x 10~^, the magnitude of u increases but in the opposite direction 
with respect to the direction of the initial spin due mainly to the retarding magnetic force. 
After t = 0.8 x 10~^, since the retarding force is switched off, the magnitude of lu starts to 
decrease by the effect of the frictional force only. For the latter case, we find from Fig. [T3] 
that the application of Hy for 10^^ s does not provide sufficient angular momentum for 
the nano top to lie down in a stable configuration. The variation of the corresponding w 
for this case is presented in Fig. [T5l 






(21) 




(22) 
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time(s) xio"' 
Figure 11: Asymptotic expansion for Uy (|22p 9 in comparison with numerical result 
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Figure 14: Nutation angle 6 for nano top when = and Hy = ifH(io - t) where = 10~^ s 
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Figure 15: Angular frequency w for nano top when = and Hy = Hli{to — t) where to = 1Q ^ s 
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